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Double-Pulse  Laser-Induced  Breakdown  Spectroscopy  of  iron  using  both  Nd:YAG  and  TEA-C02  lasers  has 
been  investigated  to  better  understand  mechanisms  of  signal  enhancement.  The  signal  dependence  on  the 
delay  between  the  two  laser  pulses  shows  an  enhanced  signal  when  the  C02  laser  pulse  interacts  with  the 
sample  before  the  Nd:YAG  pulse.  Signal  kinetics  and  a  simple  model  of  sample  heating  by  the  C02  pulse  show 
that  the  enhancement  during  the  first  700  ns  is  due  primarily  to  sample  heating.  Images  of  the  sample 
surface  after  ablation  as  well  as  time-integrated  pictures  of  the  plasma  suggest  that  particles  are  ejected  from 
the  surface  during  the  first  microseconds  after  the  arrival  of  the  C02  pulse  and  provide  fuel  for  the 
subsequent  plasma  created  by  the  Nd:YAG  laser. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Optical  emission  spectroscopy  of  laser-induced  plasmas  is  a 
powerful  spectroscopic  diagnostic  technique  for  sensing  and  other 
applications  because  of  its  ability  to  identify  all  types  of  materials  (solid, 
liquid,  gas),  at  a  close  distance  as  well  as  in  a  stand-off  configuration. 
Despite  the  increasing  popularity  of  this  technique  (also  called  LIBS  — 
Laser-Induced  Breakdown  Spectroscopy)  its  sensitivity  and  precision 
are  relatively  poor  compared  to  other  well  established  analytical 
techniques.  This  is  because  of  the  uncertainty  in  laser  energy  coupling 
to  the  sample,  significant  matrix  effects  and  relatively  high  background 
signal  in  the  atomic  and  ionic  spectra,  among  other  reasons.  One 
approach  to  improving  sensitivity  (by  increasing  the  signal-to-contin- 
uum-background  ratio)  is  the  use  of  a  double-pulse  configuration  (DP- 
LIBS).  The  aim  of  DP-LIBS  is  to  improve  the  coupling  of  the  laser  energy 
to  the  target  and  to  the  ablated  material,  leading  to  a  higher  number  of 
emitters  from  the  analyzed  sample  in  the  plasma. 

Since  the  first  demonstration  of  double-pulse  LIBS  in  1969  by 
Piepmeier  and  Malmstadt  [1],  it  has  been  demonstrated  that  two  laser 
pulses  for  LIBS  lead  to  enhanced  emission  intensities,  longer  plasma 
lifetimes  and  higher  plasma  temperatures.  To  optimize  these  effects, 
several  combinations  have  been  proposed  involving  various  geometries 
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(angle  between  the  laser  beams),  different  laser  wavelengths  and/or 
pulse  durations  and  interpulse  delay,  as  well  as  different  pulse  energies. 

A  collinear  approach  combines  two  pulses  following  the  same  path 
and  focused  upon  the  same  point  of  the  sample.  The  orthogonal 
configuration  uses  multiple  optics  to  focus  at  different  locations  (double¬ 
focus  collinear  configuration  has  also  been  demonstrated  [2]).  Two 
configurations  are  possible  with  an  orthogonal  geometry:  a  pre-ablation 
spark  to  produce  plasma  above  the  sample  prior  to  the  ablation  or  a 
reheating  spark  focused  upon  the  expanding  plasma  created  by  the  first 
pulse. 

From  femtosecond  to  microsecond  pulses,  several  combinations 
have  been  studied:  fs-fs  [3],  fs-ns  [4],  ns-ns  [5],  ns-ps  [6].  Other 
works  have  studied  the  influence  of  the  combination  of  different 
wavelengths  from  ultraviolet  (UV)  to  near-  and  mid-infrared  (NIR  and 
MIR):  UV-NIR  [5],  Visible-NIR  [8],  NIR-NIR  [5],  NIR-MIR  [6],  NIR-MW 
[9]  and  the  combination  of  pulses  with  different  energy. 

In  Babushok  et  al.  [10]  and  Scaffidi  et  al.  [4],  overviews  of  the 
different  mechanisms  leading  to  emission  enhancement  in  DP-LIBS 
experiments  are  given.  Depending  on  the  delay  time  between  the 
pulses,  several  explanations  to  the  enhancement  are  given.  From  zero- 
delay  to  tens  of  ms,  the  physics  behind  the  phenomenon  changes.  In 
Scaffidi  et  al.,  three  mechanisms  are  mainly  detailed:  Pulse-plasma 
coupling  (re-exciting  or  exciting  of  the  plasma  components),  sample 
heating  (heating  by  the  hot  plasma  close  to  the  surface  resulting  in 
increased  ablation)  and  atmospheric  effects  (creation  of  a  hot  region 
with  low  density  above  the  sample).  In  Babushok  et  al.,  an  extensive  list 
of  possible  effects  is  given  from  the  preparation  of  the  sample  by  the  first 
pulse  (smoothing,  first  crater,  onset  of  melt)  to  the  energy  coupling  with 
the  plasma  (ion  yield,  plasma  re-heating,  electron  density). 
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In  this  paper,  we  present  the  results  of  a  double-pulse  LIBS 
experiment  where  a  nanosecond  NIR  pulse  is  combined  with  a  MIR 
C02  laser  pulse  and  focused  on  an  iron  sample.  An  8-fold  enhancement 
is  reached  by  irradiating  the  sample  first  with  the  C02  laser  pulse  and 
then  ablating  it  with  the  Nd:YAG  laser.  The  effect  of  interpulse  delay 
allows  us  to  advance  a  combination  of  mechanisms  leading  to  the 
enhancement  of  the  LIBS  signal  compared  to  that  obtained  with  a 
single  Nd:YAG  laser  pulse. 

2.  Experimental  setup 

The  experimental  configuration  is  shown  in  Fig.  1.  Two  lasers  are 
fired  onto  the  same  spot  of  a  sample  and  the  light  emitted  by  the 
plasma  is  collected  and  guided  to  a  spectrometer.  The  spatial  overlap 
of  the  two  pulses  is  shown  on  the  inset  in  Fig.  1. 

2.1.  Laser  sources 


symmetric  resonator  and  emits  at  10.6  pm.  An  intra-cavity  iris  is  used 
to  limit  the  number  of  higher  order  transverse  modes  and  provide  a 
circularly  symmetric  beam  profile  as  viewed  using  an  IR  pyroelectric 
beam  camera  (Spiricon  Pyrocam  III).  The  laser  energy  is  continuously 
tunable  between  70  pj  and  80  mj  by  using  a  set  of  crossed  polarizers. 
The  pulse  energy  is  measured  using  the  same  pyroelectric  joulemeter 
described  above.  The  peak-to-peak  and  RMS  fluctuations  in  pulse 
energy  are  3.5%  and  0.98%,  respectively.  The  temporal  profile  is 
measured  with  a  multiple  junction  photovoltaic  detector  (Vigo 
Systems  PVM-10.6)  with  nanosecond  response  time.  The  temporal 
profile  is  composed  of  a  100  ns  pulse  followed  by  a  1  ps  tail  (Fig.  7a). 
The  laser  is  focused  with  a  26  cm  focal  length  lens  at  13  degrees  from 
normal,  and  the  lens-to-sample  distance  is  25.5  cm.  The  beam  radius 
at  the  sample  location  is  approximately  280  pm  as  determined  using  a 
knife-edge  scan.  An  irradiance  of  30xl06  W/cm2  (assuming  a  pulse 
duration  of  1  ps)  was  obtained  with  the  laser  operating  at  75  mj  per 
pulse. 


The  Nd:YAG  laser  (Quantel  Brilliant)  is  Q-switched  and  operated  at 
the  fundamental  wavelength  of  1064  nm.  It  is  externally  triggered  at 
0.33  Hz  as  limited  by  the  C02  laser.  A  thin-film  polarizer  is  used  in 
conjunction  with  a  \/2  plate  to  provide  continuous  control  of  the  laser 
energy  from  5  mj  up  to  340  mj  without  changing  the  flashlamp 
voltage  or  the  Q-switch  delay  time.  The  pulse  energy  is  measured  with 
a  calibrated  pyroelectric  joulemeter  (Gentec  SOLO  2  with  a  QE-25 
head).  The  peak-to-peak  and  RMS  fluctuations  in  pulse  energy  are 
1.51%  and  0.44%  respectively.  The  full  width  half  maximum  (FWHM) 
pulse  duration  is  5  ns  as  measured  with  a  fast  photodiode  (Thorlabs) 
placed  behind  a  dielectric  turning  mirror.  The  laser  is  focused  with  a 
13.5  cm  focal  length  lens  at  normal  incidence,  and  the  lens-to-sample 
distance  is  13  cm.  The  beam  radius  at  the  sample  location,  determined 
by  the  knife-edge  scan  technique,  is  approximately  230  pm.  The  Nd: 
YAG  laser  was  operated  at  39  mj  per  pulse,  corresponding  to  an 
irradiance  of  4.7  x  109  W/cm2. 

The  C02  laser  is  a  transversely  excited  atmospheric  (TEA)  C02  laser 
with  components  made  by  Lumonics.  It  utilizes  an  external  1  m  half- 


A/2  Plate 


Fig.  1.  Experimental  setup. 


2.2.  Collection  optics  and  spectrometer 

The  plasma  emission  is  collected  using  an  f/1  quartz  lens  and  is 
coupled  to  a  1  m  round-to-line  UV-grade  fiber  bundle  using  a  second 
(f/2)  lens.  The  collection  optics  are  at  40°  from  the  surface  normal.  The 
fiber  bundle  is  f-matched  to  a  0.5  m  Czerny-Turner  imaging 
spectrometer  (Acton  2500i,  Princeton  Instruments)  with  a  1800 
lines/mm  grating.  The  light  is  detected  using  an  intensified  charge 
coupled  device  (Andor  iStar  DH720-25F-03)  with  1024x256  pixels 
and  a  25  mm  Gen  II  intensifier.  The  intensifier  delay  is  measured  from 
the  rising  edge  of  the  second  pulse.  Each  spectrum  is  the  average  of  30 
single  acquisitions. 

2.3.  Timing 

The  timing  of  the  experiment  is  controlled  using  two  digital  delay 
generators  (Stanford  Research  Models  DG-645  and  DG-535).  The  DG- 
645  triggers  the  flashlamp  of  the  Nd:YAG  laser  and  also  triggers  the 
DG-535  at  10  Hz.  The  DG-535  divides  the  frequency  to  0.33  Hz  and 
triggers  the  Q-switch  of  the  Nd:YAG  laser,  the  C02  laser  and  the  ICCD 
camera.  Timing  uncertainty  between  Nd:YAG  and  C02  pulses  is 
primarily  limited  by  the  jitter  in  the  C02  laser  pulse  onset  and  has 
peak-to-peak  fluctuations  of  approximately  42  ns  and  a  standard 
deviation  of  9  ns.  The  timing  configuration  for  the  interpulse  delay 
and  the  acquisition  is  shown  in  Fig.  2. 

2.4.  Samples 

The  target  samples  were  a  100  pm-thick  foil  of  99.9  +  %  pure  iron 
(Product  #:  356808,  Sigma-Aldrich,  St.  Louis,  MO,  USA).  Each  foil  was 
cleaned  with  acetone  and  mounted  so  that  it  remained  flat  and 
vertical  to  minimize  any  deposition  of  material  after  the  ablation.  The 


Fig.  2.  Timing  configuration  of  the  experiment.  The  interpulse  delay  Tp  is  set  as  positive 
when  the  Nd:YAG  laser  precedes  the  C02  laser  pulse  and  negative  in  the  opposite  order. 
The  acquisition  delay  Td  is  set  to  1  ps  after  the  second  pulse  and  the  duration  of 
acquisition  Tw  to  25  ns. 
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sample  was  translated  between  each  double  pulse  laser  irradiation  in 
order  to  present  a  fresh  location  for  each  shot. 

3.  Results  and  discussion 

The  original  intention  of  this  double-pulse  LIBS  study  was  to  use 
10.6  pm  radiation  from  the  C02  laser  to  re-heat  or  raise  the 
temperature  of  the  plasma  created  by  the  Nd:YAG  laser  after  it  had 
diffused  away  from  the  surface.  Since  the  critical  density  of  a  plasma, 
the  point  at  which  light  is  reflected,  scales  as  the  square  of  the 
wavelength, 

cr  _  f2nc\2e0me 
e  \  e  )  \2 

longer  wavelength  light  will  interact  with  the  plasma  farther  down 
the  density  gradient  in  the  corona  of  the  plasma  produced  by  the  Nd: 
YAG  laser.  With  the  relatively  long  scale  lengths  produced  in  the 
plasma,  inverse  Bremsstrahlung  absorption  (IBA)  occurs  at  densities 
well  below  the  critical  density. 

Previous  studies  using  multiple  wavelengths  for  double-  pulse  LIBS 
has  been  reported  [5]  using  both  UV  and  NIR  pulses,  and  a  correlation 
was  found  between  signal  enhancement  and  plasma  temperature  [5] 
showing  that  the  signal  increases  with  the  temperature. 

A  comparison  between  single-pulse  and  double-pulse  LIBS  spectra 
is  shown  in  Fig.  3.  Each  spectrum  represents  the  average  of  30 
separate  shots  obtained  with  1  ps  delay  (measured  with  respect  to  the 
second  laser  pulse)  and  an  acquisition  time  of  25  ns.  This  acquisition 
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time  was  chosen  for  determining  instantaneous  plasma  temperature. 
Spectra  for  the  single-pulse  case  of  Nd:YAG  and  C02  are  shown  in  a) 
and  b)  respectively,  and  the  double-pulse  case  for  Nd:YAG  preceding 
the  C02  pulse  by  1  ps  and  Nd:YAG  following  the  C02  pulse  by  1  ps  are 
shown  in  c)  and  d)  respectively.  In  all  cases,  the  energy  of  Nd:YAG  and 
C02  lasers  were  23  J  cm-2  and  30  J  cm-2  respectively. 

The  spectral  lines  in  Fig.  3a  are  of  singly-ionized  iron.  The  weak 
signal  in  Fig.  3b  can  be  attributed  to  minimal  ablated  mass  and  low 
plasma  temperature  during  the  relatively  long  interaction  time  of  the 
C02  pulse.  For  the  double  pulse  case,  with  the  C02  laser  pulse 
reheating  the  plasma  (Fig.  3c),  there  is  no  signal  enhancement  relative 
to  the  signal  acquired  from  the  Nd:YAG  only.  This  can  be  partially 
attributed  to  the  acquisition  delay  time  of  1  ps  measured  from  the 
second  laser  pulse;  therefore,  2  ps  after  the  Nd:YAG  pulse.  This 
measurement  does  not  imply  significant  plasma  heating  by  the  C02 
pulse.  However,  the  effects  of  plasma  reheating  using  a  C02  laser  pulse 
(increased  line  intensities)  have  been  demonstrated  previously. 

In  work  by  Killinger  et  al.  [6],  an  enhancement  of  both  atomic  (up 
to  a  factor  of  60)  and  ionic  (up  to  a  factor  of  300)  aluminum  lines  for  a 
high  purity  alumina  sample  was  observed  when  a  Nd:YAG  pulse  (5  ns, 
50  mj  mm-2)  was  followed  by  a  C02  pulse  (100  ns  peak  and  5  ps  tail, 
70  mj  mm-  2).  The  difference  in  behavior  compared  to  our  results  can 
be  attributed  to  both  the  different  energy  densities  used,  since  the 
signal  from  a  nanosecond  laser  plasma  is  approximately  proportional 
to  I3/2  [7],  I  being  the  laser  intensity,  and  sample  properties.  A  metal 
sample  was  used  in  our  case  and  an  insulator  in  their  case,  each  having 
different  thermodynamic  properties.  A  complete  comparison  cannot 
be  made  with  our  work  since  the  timing  configuration  is  not 
mentioned  in  their  publication. 

To  understand  the  effect  of  the  C02  pulse  in  double-pulse  LIBS,  the 
time  delay  between  pulses  has  been  varied  from  the  C02  pulse 
preceding  the  Nd:YAG  pulse  by  5  ms  to  the  C02  pulse  following  Nd: 
YAG  pulse  by  10  ps.  For  clarity,  Fig.  4  shows  the  intensity  of  the  Fe(II) 
274.65  nm  ionic  line  for  a  delay  between  —  lOps  and  +  10  jus.  A 
negative  delay  means  the  C02  pulse  reaches  the  sample  before  the  Nd: 
YAG  pulse.  Although  there  is  no  plasma  reheating  by  the  C02  pulse,  a 
signal  enhancement  of  approximately  8  times  occurs  for  an  interpulse 
delay  of  —  1  ps.  This  signal  enhancement  is  closely  correlated  to  an 
enhancement  in  plasma  temperature,  as  calculated  using  a  Boltzmann 
plot  [11]  with  the  parameters  listed  in  Table  1.  The  temperature  for  a 


Interpulse  delay  (ps) 


Fig.  3.  LIBS  spectra  of  iron  from  a)  a  single  Nd:YAG  pulse,  b)  a  single  C02  laser  pulse, 
c)  double  pulse  spectra  with  C02  following  Nd:YAG  by  1  ps  compared  with  Nd:YAG 
alone  LIBS  spectrum  and  d)  double  pulse  spectra  with  C02  preceding  Nd:YAG  by  1  ps. 
(Acquisition  delay:  1  ps,  duration  of  acquisition:  25  ns).  In  a),  the  dotted  line  shows  a 
background  spline  interpolation  and  the  stars  indicate  the  lines  used  for  plasma 
temperature  calculation.  Baseline  correction  is  applied  for  all  spectra. 


Fig.  4.  Background-corrected  peak  intensity  of  the  Fe  II  274.65  nm  line  as  a  function  of 
the  interpulse  delay.  A  negative  delay  means  the  C02  pulse  reaches  the  samples  before 
the  Nd:YAG  pulse.  A  positive  delay  is  for  the  opposite  timing  configuration.  Comparison 
is  made  with  the  single-pulse  cases.  The  error  bars  represent  the  standard  deviation 
over  30  single  acquisitions.  The  inset  in  the  upper  right  corner  is  a  zoom  between  —  3  ps 
and  0  ps  showing  the  maximum  enhancement  at  —  1  ps. 
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Table  1 

Parameters  of  the  spectral  lines  used  in  the  plasma  temperature  calculations  [12]. 


Wavelength  \ui 
(nm) 

Upper  level  energy  Eu 
(eV) 

Degeneracy 

gu 

Transition  probability  Au[ 

(s'1) 

268.48 

8.43 

10 

1.4  108 

269.26 

8.37 

12 

1.2  108 

270.40 

7.97 

8 

1.2  108 

271.44 

5.55 

6 

5.5  107 

273.96 

5.51 

8 

1.9  108 

274.65 

5.59 

6 

1.9  108 

274.93 

5.55 

8 

2.1  108 

275.57 

5.49 

10 

2.1  108 

278.37 

7.70 

10 

7.0  107 

single-pulse  Nd:YAG  plasma  is  15  000  ±2000  K  and  23  000  ±2000  K 
for  double  pulse  with  —  1  ps  interpulse  delay. 

When  the  C02  pulse  is  fired  before  the  Nd:YAG,  two  mechanisms 
may  occur:  (i)  the  C02  pulse  heats  the  sample  before  the  arrival  of  the 
Nd:YAG  pulse,  (ii)  the  Nd:YAG  interacts  with  the  weak  plasma  created 
by  the  C02  pulse,  ablating  the  sample  softly  and  extracting  non 
excited  particles  from  the  sample  with  size  distribution  ranging  from 
atomic  to  clusters. 

3.1.  Heating  of  the  sample  by  the  C02  pulse 

Several  studies  have  shown  a  positive  correlation  between  the 
intensity  of  the  LIBS  signal  intensity  and  the  sample  temperature 
[13,14].  This  phenomenon  is  explained  by  less  energy  being  required 
to  melt  or  vaporize  the  sample  —  the  important  parameters  in  long- 
pulse  ablation  above  room  temperature.  Furthermore,  the  reflectivity 
of  metals  decreases  when  the  surface  temperature  increases  [15]. 
These  two  factors  could  lead  to  a  better  energy  coupling  between  the 
Nd:YAG  pulse  and  the  sample  if  previously  heated  by  the  C02  pulse.  To 
understand  the  importance  of  this  phenomenon,  a  simple  model  of 
heat  conduction  is  used. 

The  model  is  one-dimensional  (ID)  since  the  absorption  length  of 
a  radiation  with  a  10.6  pm-wavelength  at  298  K  in  the  target  is 
approximately  30  nm  (la  =  l/a,  with  a  =  3.3  107m-1),  which  is 
much  smaller  than  either  laser  beam  diameter.  The  ID  heat  equation 
has  the  following  form: 


dT(x ,  t) 
~~dF~ 


d  I"  /  k  \  dT (x,  t) 
dx  yCpp J  dx 


+<5sm 


given  in  Table  2.  The  temporal  profile  of  the  heat  source  I0(t)  is  an 
experimental  profile  acquired  from  a  75  mj  C02  pulse.  The  model 
yields  an  estimate  of  the  temperature  on  the  sample  surface  and 
doesn't  include  plasma  shielding  or  ablation  of  the  sample  by  the  laser 
(the  energy  at  the  surface  is  considered  completely  transformed  into 
heat).  The  lack  of  absorption  of  the  laser  pulse  by  the  plasma  leads  us  to 
consider  that  this  simple  model  will  give  a  higher  estimate  of  the 
temperature  at  the  surface  of  the  sample.  However,  no  conclusions  on 
the  ablation  products  or  the  plasma  characteristics  can  be  extracted. 

The  calculated  temperature  distribution  in  the  target  is  plotted  as  a 
function  of  time  as  shown  in  Fig.  5.  The  first  conclusion  is  that  the 
melting  point  Tm  ( 181 1  K  for  iron)  is  not  reached  with  the  irradiance  of 
the  C02  laser  pulse  used.  The  maximum  temperature  at  the  surface  is 
approximately  860  K  and  occurs  at  550  ns.  This  is  supported  by  the 
microscope  images  of  the  sample  irradiated  by  only  the  C02  laser 
(Fig.  6b)  where  small  degradation  of  the  material  by  the  C02  pulse  is 
observed  after  30  shots  on  the  same  spot.  The  crater  depth  after  30  shots 
(using  a  Zygo  NewView  6300  3D  optical  profiler)  is  not  discernible 
within  the  roughness  of  the  sample  by  the  instrument  and  is  smaller 
than  the  absorption  length  (30  nm)  which  is  indicative  of  nanometric 
layers  being  ablated  by  each  C02  pulse.  Given  minimal  material  removal 
by  melting  and  vaporization,  the  Nd:YAG  pulse  following  the  C02 
ablation  is  not  influenced  by  significant  sample  deformation. 

Given  the  maximum  temperature  occurs  at  t=  550  ns,  a  maximum 
signal  enhancement  at  an  interpulse  delay  of  1  ps  is  not  completely 
explained  by  sample  temperature  effects. 

3.2.  Interaction  of  the  Nd:YAG  laser  pulse  with  C02  laser  induced  plasma 

With  little  apparent  sample  degradation  as  a  result  of  the  C02  laser, 
a  continuum-like  emission  is  detected  as  shown  in  Fig.  3b. 

This  continuum  can  be  attributed  to  bremsstrahlung  and/or  poly¬ 
atomic  clusters.  These  emitters  are  extracted  by  the  C02  laser  pulse  during 
the  interaction  with  the  surface,  when  a  nanometric  layer  is  ablated. 
When  the  laser  removes  material  from  the  surface,  the  size  distribution  of 
particles  is  limited  by  the  absorption  depth  1/a  [1 7].  Two  types  of  emitters 
are  localized  in  time-integrated  pictures  of  the  plasma  in  the  visible  region 
shown  in  Fig.  7b  (using  a  manually  triggered  Nikon  D40x  DSLR  with 
300  mm  1:1  macro  zoom  lens  at  f/8  with  1  to  2  second  exposure  time  to 
integrate  the  entire  plasma  lifetime).  Light  emission  is  observed  at  two 
different  locations:  one  component  along  the  normal  and  the  other  close 
to  the  surface.  The  paraxial  emission  can  be  attributed  to  small  emitters 
(atoms  and  ions),  ejected  with  high  speed  from  the  sample  (tens  of 


with  S(x,t)  =  /0(t)e_ax(l  —R)-Io(t)  is  the  laser  irradiance  (in  Won-2) 
and  the  other  parameters  are  listed  in  Table  2.  This  calculation  is  done 
via  a  matricial  implicit  algorithm  for  solving  parabolic  partial 
differential  equations  [16].  The  conditions  of  the  algorithm  are  also 


Table  2 

Physical  and  numerical  parameters  used  in  our  ID  heat  equation  model. 


Parameters 

Values 

Thermal  conductivity,  k  (W  cm-1  K-1) 

0.802 

Specific  heat,  CP(J  g-1  K-1) 

0.449 

Mass  density,  p(g  cm-3) 

7.87 

Absorption  coefficient,  a  (m-1)  (\=  10.6  pm) 

3.3  107 

Reflectivity,  R  (X  =  10.6  pm) 

0.967 

Melting  point,  Tm  (K) 

1811 

Space  range 

[x0  =  0  pm;  Xf=  1 0  pm] 

Space  sampling 

701  points 

Time  range 

[t0=  —650  ns;  tf=  1849  ns] 

Time  sampling 

2500  points 

Initial  condition 

T(x,t=0)  =  300/C 

Boundary  conditions  on  x0  and  Xf 

dT(x  =  Xo,t)  _  q  dT(x  =  xf,t)  _  Q 
dx  ’  dx 

The  value  corresponding  to  the  solid  form  is  given  by  the  literature  [19]. 


Fig.  5.  Calculated  temperature  distribution  in  the  target  as  a  function  of  time  for  an  iron 
sample  after  the  irradiation  with  the  laser  pulse  depicted  in  Fig.  7a.  (see  physical  and 
algorithmic  parameters  in  Table  1). 
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Fig.  6.  Microscope  pictures  of  the  iron  sample  after  a)  a  single  Nd:YAG  laser  shot,  b)  an  accumulation  of  30  C02  laser  shots,  c)  a  combination  of  a  Nd:YAG  laser  shot  followed  by  a  C02 
laser  shot  1  ps  later  and  d)  a  combination  of  a  C02  laser  shot  followed  by  a  Nd:YAG  laser  shot  1  ps  later. 


km  s“ 1 )  [18],  whereas  the  marginal  emission  can  be  due  to  slow  heavy 
emitters  (clusters)  (hundreds  of  m  s-1)  [18].  This  spatially  divided 
emission  is  compared  with  the  light  emission  of  the  plasma  created  by 
the  Nd:YAG  laser  only  (acquired  in  the  same  conditions)  where  the  plasma 
is  homogeneously  confined  to  the  surface  and  its  expansion  is  comparable 
to  the  beam  waist. 

The  spatially  integrated  temporal  trace  of  the  visible  emission  is 
shown  in  Fig.  8c.  The  first  peak  is  due  to  bremsstrahlung  from  electronic 
collisions  with  fast  emitters,  and  is  shown  by  the  paraxial  emission  in 
Fig.  7b.  During  the  course  of  the  C02  laser  pulse,  there  is  a  second 
emission  increase  peaking  at  800  ns.  This  is  attributed  to  heavy  particles 
extracted  from  the  surface  during  the  relatively  long  C02  pulse 
evolution.  The  low  velocity  of  these  emitters  restricts  their  position  to 
within  hundreds  of  microns  of  the  sample,  as  seen  in  Fig.  7b. 

When  the  Nd:YAG  pulse  arrives  1  ps  after  the  C02  pulse,  a 
maximum  coupling  with  the  C02  induced  plasma  occurs  because  the 
maximum  number  of  particles  has  been  extracted  from  the  sample 
(since  there  is  no  excitation  by  the  C02  laser  anymore).  Particle 
ablation  by  the  Nd:YAG  laser  is  efficient  because  of  their  high  surface- 
to-mass  ratio,  and  at  this  time  an  enhanced  ionic  LIBS  spectrum  is 
observed  (Fig.  3d)  compare  to  the  surface  ablation  by  the  Nd:YAG  only 
(Fig.  3a).  There  is  a  significant  difference  in  the  visible  emission  for  the 
case  of  C02  pulse  preceding  the  Nd:YAG  pulse  by  1  ps  and  the  C02 
pulse  following  the  Nd:YAG  pulse  by  1  ps  as  reflected  in  Fig.  7c  and  d. 
These  images  show  a  less  intense  surface  plasma  in  the  second  case 
suggesting  the  absorption  of  the  Nd:YAG  above  the  surface.  This  is 
also  observed  in  the  crater  images  shown  in  Fig.  6.  The  edges  of  the 
crater  induced  by  the  Nd:YAG  in  the  double  pulse  case  of  C02  first  are 
smoother  than  in  the  case  when  a  single  Nd:YAG  pulse  reaches  the 
sample  first  (Fig.  6a  and  c). 


3.3.  Combination  of  mechanisms  for  LIBS  signal  enhancement 

Maximum  signal  enhancement  occurs  when  the  C02  precedes  the 
Nd:YAG  pulse  by  1  ps.  This  condition  is  shown  in  Fig.  8  along  with  the 
normalized  plasma  emission  and  calculated  surface  temperature.  We 
can  divide  the  timescale  in  two  stages: 

3.3.1.  The  first  700  ns 

During  the  first  700  ns,  there  is  a  strong  correlation  between  the 
surface  temperature  and  the  plasma  emission.  The  100  ns  peak  of  the 
C02  laser  pulse  quickly  heats  the  sample  from  approximately  300  K  to 
650  K.  The  tail  of  the  pulse  reaches  a  hot  surface  and  heats  the  sample 
to  860  K.  The  signal  follows  the  same  trend  until  approximately 
700  ns.  This  enhances  the  production  of  emitting  products.  These 
products  are  mainly  heavy  emitters  rather  than  ions  or  atoms 
common  to  most  LIBS  plasmas  and  there  is  little  surface  modification, 
leading  to  the  conclusion  that  these  heavy  particles  have  nanometer 
size.  When  the  Nd:YAG  pulse  reaches  a  hot  surface,  less  energy  is 
needed  for  ablation  resulting  in  increased  mass  ablation  and  LIBS 
signal. 

3.3.2.  Between  700  ns  and  1 0  ps 

The  extraction  of  particles  from  the  surface  increases  until 
approximately  0.7  to  1  ps.  At  this  time,  the  coupling  between  the 
C02  induced  plasma  and  the  Nd:YAG  pulse  is  maximum.  Emission 
from  the  ablation  of  particles  with  a  distribution  of  velocity  around 
500  m  per  second  will  remain  within  the  field  of  view  (approximately 
2  mm)  for  approximately  4  ps.  However,  noticeable  enhancement  of 
the  LIBS  signal  occurs  until  10  ps.  Further  experiments  are  in  progress 
to  better  explain  the  mechanisms  of  such  a  long  enhancement. 
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Fig.  7.  Time-integrated  pictures  of  the  plasma  produced  by  a)  a  single  Nd:YAG  laser  shot,  b)  a  single  C02  laser  shot,  c)  a  combination  of  an  Nd:YAG  laser  shot  followed  by  a  C02  laser 
shot  1  ps  later  and  d)  a  combination  of  a  C02  laser  shot  followed  by  a  Nd:YAG  laser  shot  1  ps  later.  The  arrows  represent  the  direction  of  the  incoming  laser  and  the  dotted  lines  the 
laser  beam  diameter  2w0. 


4.  Conclusion 

The  optical  emission  from  laser-induced  plasma  in  a  double-pulse 
LIBS  experiment  (Nd:YAG  at  1064  nm  and  C02  laser  at  10.6  pm)  on  iron 
has  been  studied  to  understand  which  mechanisms  are  responsible  for 
signal  enhancement  in  this  configuration.  The  C02  laser  pulse  heats  the 
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Fig.  8.  Optimal  timing  for  the  maximum  of  LIBS  signal  enhancement:  a)  C02  laser  pulse 
intensity,  b)  Nd:YAG  laser  intensity,  c)  the  plasma  emission  from  the  C02  pulse  only  and 
d)  the  calculated  surface  temperature  of  the  sample. 


sample  before  the  arrival  of  the  Nd:YAG  laser  pulse  during  the  first 
700  ns,  without  melting  the  sample  or  creating  noticeable  surface 
deformation.  Nevertheless,  the  bright  continuum  emission  representa¬ 
tive  of  heavy  and  slow  moving  particles  extracted  and  ejected  during  the 
C02  laser  interaction  with  the  iron  sample.  These  particles  act  as  fuel  for 
the  second  pulse,  and  lead  to  greatest  signal  enhancement  approxi¬ 
mately  one  microsecond  after  the  beginning  of  the  C02  laser  pulse. 
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